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Abstract: The determination of heavy element abundances from planetary nebula (PN) spectra pro- 
vides an exciting opportunity to study the nucleosynthesis occurring in the progenitor asymptotic giant 
branch (AGB) star. We perform post-processing calculations on AGB models of a large range of mass 
and metallicity to obtain predictions for the production of neutron-capture elements up to the first 
s-process peak at strontium. We find that solar metallicity intermediate-mass AGB models provide a 
reasonable match to the heavy element composition of Type I PNe. Likewise, many of the Se and Kr 
enriched PNe are well fitted by lower mass models with solar or close-to-solar metallicities. However 
the most Kr-enriched objects, and the PN with sub-solar Se/O ratios are difficult to explain with AGB 
nucleosynthesis models. Furthermore, we compute s-process abundance predictions for low-mass AGB 
models of very low metallicity (Fe/H w —2.3) using both scaled solar and an a-enhanced initial com- 
position. For these models, O is dredged to the surface, which means that abundance ratios measured 
relative to this element (e.g., X/O) do not provide a reliable measure of initial abundance ratios, or of 
production within the star owing to internal nucleosynthesis. 

Keywords: stars: AGB and post-AGB — abundances — ISM: abundances — planetary nebulae: 
general 



1 Introduction 



After the thermally-pulsing AGB (TP- AGB) phase is 
terminated, low to intermediate-mass stars (~ 0.8 to 
8Mq) evolve at near constant luminosity to become 
post-AGB objects. If the ejected envelope has suffi- 
cient time to become ionized by the hot central star 
before dissipating, then the object will also be ob- 
served as a planetary nebula (PN). For recent r eviews 
of TP-AGB and post- AGB stars, see iHerwid (|2005h 
and lvan W inckcl ( 20Q3) , respectively. The ionised neb- 
ula is comprised of material from the convective en- 
velope that once surrounded the core, hence, nebu- 
lar abundances may reveal information about the effi- 
ciency of mixing events and chemical processing that 
took place during previous evolutionary phases, in ad- 
dition to the initial composition of the p arent star (e.g., 
iDopita et all I 19971 : iMaerini et al.ll2009h . 

During the life of a star, mixing between the nuclear- 
processed core and envelope may occur, and this has 
the effect of enriching the surface composition in the 
products of H-burning (first and second dredge-up; hot 
bottom burning), and He- burning (the third dredge- 
up). Note that the first dredge- up occurs after core 
H-exhaustion, whereas the second dredge-up after core 
He-exhaustion. Both the third dredge up (TDU) and 
hot bottom burning (HBB) occur during the TP-AGB 
phase . We refer the reader to, e.g.. Hoothrovd fc Sackmannl 
(| 19991 ) for details on the first and second dredge up. 

A TP-AGB star is characterised by two nuclear 
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Figure 1: Schematic structure of an AGB star. 
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burning shells above a degenerate C-O core, surrounded 
by a deep convective envelope (see Fig. [TJ. The TDU 
mixes material from the He-intershell to the envelope. 
This material is composed primarily of 4 He and 12 C, 
along with trace quantities of heavy elements produced 
by the slow neutron capture process (the s-process). 
The TDU may occur after each instability of the he- 
lium shell, slowly enriching the envelope in carbon 
and heavy elements. HBB occurs in intermediate- mass 
AGB stars (~ 4 to 8M0) when the base of the con- 
vective envelope can dip into the top of the hydrogen 
burning shell. 

Abundances derived from PN spectra provide a 
complimentary data set to the abundances derived from 
the spectra of cool evolved stars. For example, the 
elemental abundances of He, Ne, and Ar can be ob- 
tained, along with the abundances of C, N, O, S, and CI 
(lAller fc Czvzaklll983l ; IStanghellini et alll2000l ; iLeisv fc Denrjep. 
l2006h . Recent observations have also revealed that 




some PNe are enriched in heavy elements that can 
be produce d by the .s-process including Ge, Se, Kr, 
Xe, a nd Ba (jSharpee et al ]|2007l : ISterling fc Dinersteinl 
|2008| . hereafter SD08). SD08 obtained Se and Kr abun- 
dances for 120 Galactic PNe, and investigated trends 
between s-process enrichments and PN morphology 
and other nebular and stellar characteristics. It was 
found that while some PN have large enrichments of 
Se and Kr (e.g., [Se,Kr/0] > 1), Type I PN have lower 
s-process abundances, on average, than the sample as 

a whole. 

In IKarakas et all ( 2009) we compared nucleosyn- 
thesis predictions from models of intermediate-mass 
AGB stars to the results from SD08. Here we briefly 
summarize the results from Karakas et al. (2009), and 
provide new results for lower mass AGB models of so- 
lar and halo compositions. 



2 Type I Planetary Nebula 



The aims of IKarakas et all (2009) were to compare 
the s-process predictions from intermediate-mass AGB 
stars to the Se and Kr compositions of the Type I 
PN identified in the SD08 sample. Such a comparison 
may be used to constrain the efficiency of the TDU 
in AGB models, which is not well determined. Also 
this information may help elucidate to what extent 
intermediate-mass AGB stars contribute to the Galac- 
tic inventory of s-process elements. The motivation to 
use Type I PN comes from their high He/H and N/O 
ratios characteristic of p roton-capture nucleosynthesis 
( Stanghcll ini et al"1 l2006l as well as spatial distribu- 
tions and kinematics indi cative of a young population 
(jCorradi fc Sch warz 1995) with initial masses between 
w 2.5 to 6M . 

Furthermore, we had also hoped to identify the 
neutron source(s) operating in intermediate-mass AGB 
stars. In low-mass stars it has been shown that the 
13 C(a,n) 16 reaction is the dominant source of neu- 
trons, where as the 22 Ne(a, n) 25 Mg reaction only plays 
a minor role l|Lambert et al.lll995l ; lGallino etaflll998ft 
The 22 Ne(a, ri) reaction has been suggested as the dom- 
inant neutron source in intermediate-mass AGB stars 
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Figure 2: The surface abundances of neutron-capture 
elements around the first s-process peak. The solid line 
indicates abundances from the 3M0, Z = 0.02 model, 
.shed line abundances from the 6.5Mq, Z = 0.02 
moaei, the dot-dashed line abundances from the 5M© , 
Z = 0.008, and the dotted line abundances from the 
5Mq, Z = 0.004 model. Abundances are taken at the 
tip of the TP- AGB, after the last computed thermal 
pulse. The boxes represent the range of observed Se 
and Kr abundances for Type I PNe, see the text for 
details. 



(jTruran fc Ibenlll977l ; iGarci'a-Hernandez et all 120061 ). 
The s-process abundance pattern that results from 
each reaction is quite different, owing to differing time- 
scales over which the neutrons are released along with 
maximum neutron densities that vary by orders of mag- 
nitude: ~ 10 8 for the 13 C reaction compared to ~ 10 13 
neutrons per cm -3 for 22 Ne. One important result, 
for example, is the prediction of low Rb/Sr ratios in 
low-mass AG B stars, which has b een observationally 
verified fe.g.. iLambert et al.lll995l ). 

The details of the stellar models and nucle osynthe- 

sis cal culat ions have been discuss e d in d etail in lKarakas et al.l 
(2009) and lKarakas fc Lattanzid (|2007l ). Here we sum- 
marize the main findings. It is important to point out 
that only models with HBB (> 4Mq) would be ob- 
served as Type I PNe, whereas observations put the 
minimum mass at ~ 2.5M©. Clearly, the lowest mass 
range would require either rapid stellar rotation, ex- 
tra mixing, and/or the effect of a binary companion to 
assist in producing the Typ e I composition. 

In lKarakas et all (|2009t ) we found that the HBB 
AGB models of solar or near solar metallicity (in par- 
ticular Z — 0.02 and 0.008) were a reasonable match 
to the Se and Kr abundance distribution of the Type I 
PNe, see Fig. [2] In Fig. [2]the abundances are shown as 
[X/Ar] ratios, and are plotted as a function of proton 
number, Z. The hatched boxes around the elements 
Se and Kr indicate the range of the Type I PN abun- 
dances from the SD08 sample, with the solid vertical 
line extending the boxes by an extra 0.3 dex, indicative 
of the typical errors. 

The intermediate-mass low metallicity models (Z = 
0.004) produced Se and Kr abundances that were in 
excess of the Type I PN abundances (see the dotted 
line in Fig. [2]), indicating that the PN we see today 
did not form from such low-metallicity objects. This 
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result is consistent with the PN evolving from a young, 
disk population. The 3M0 models of solar metallicity 
were a marginal fit to the Se and Kr data for Type 
I PNe. The 3Mq models became carbon rich and we 
would not predict them to be observed as Type I ob- 
jects with high He/H and N/O ratios. Rapid rotation 
during the main sequence could, however, bring some 
CN-processed material to the surface. 

Binarity may be the condition nec essary for shap- 
ing b ipolar PNe (and perhaps all PNe) (|Moe fc De Marcq 
l2006h . Many, but not all, Type I PNe are also bipolar 
in shape. If binary interactions can produce a Type I 
composition, then this may lead to the TP-AGB phase 
being terminated early, before many thermal pulses 
have occurred. Hence the prediction would be for lower 
s-process elements to be produced, on average, com- 
pared to star that evolve to the tip of the AGB phase. 
This would also be consistent with the observations of 
no s-process enrichment observed in Type I PNe. 

3 Se and Kr abundances in low- 
mass AGB stars 

In Fig. [3] we show the surface abundances of Zn, Ge, 
Se, Br, Kr, and Sr at the tip of the AGB phase for 
three low-mass AGB models. The solid line represents 
the abundances from the 3M0, Z = 0.02 model, the 
dashed line abundances for the 2.IM0, Z — 0.008 
model, and the dot-dashed line abundances for the 
2.5M0, Z — 0.008 model. Abundances are shown as 
[X/O] ratios. The hatched boxes around the elements 
Se and Kr indicate the range of the PN abundances 
for the full sample from SD08, where the solid vertical 
line extending the boxes by an extra 0.3 dex, indica- 
tive of the typical errors. All models include a partial 
mixing zone of O.OO2M0. This mixing zone produces a 
13 C pocket during the interpulse period, allowing the 
reaction 13 C(a,n) le O to release free neutrons. Note 
that models without partial mixing zo nes result in no 
enric hment of s-process eleme nts (e.g.. | Karakas et al.l 
200711. We refer the reader to lBusso et al.l ljl999h and 
Herwid (|2005l ) br detailed discussion about the forma- 
tion of 13 C pockets and related uncertainties. 

From Fig. [3] we see that the solar metallicity and 
Z — 0.008 models are within the observed range of Se 
and Kr PN abundances. Starting with a scaled solar 
abundance distribution, with [s/O] = 0, it is not pos- 
sible to obtain sub-solar Se abundances, i.e., [Se/O] < 
0. Likewise, the most Kr-enriched objects, with [Kr/O] 

> 1, are also not matched by the models. 

The boxes indicating the range of derived PN abun- 
dances are somewhat misleading, because they do not 
convey the Se and Kr abundance distribution of PN, 
and most importantly, the number with extreme Se 
and Kr abundances. Fig. 4 from SD08 shows that 
there are only six PN out of 120 with low Se abun- 
dances, that is [Se/0,Ar] < —0.3. There are also only 
six PN with high Kr abundances, that is [Kr/0,Ar] 

> 1. Indeed, the rest of the sample of 120 PNe shows 
[Se/0,Ar] between and 1, consistent with the stellar 
models. Likewise, most of the Kr-enriched PN have 
< [Kr/O, Ar] < 1.0, consistent with the AGB mod- 
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Figure 3: The surface abundances of neutron-capture 
elements around the first s-process peak. The solid line 
indicates abundances for the 3M©, Z = 0.02 model, 
the dashed line abundances for the 2.IM0, Z — 0.008 
model, and the dot-dashed line abundances for the 
2.5M0, Z = 0.008 model. Abundances are taken at 
the tip of the TP-AGB, after the last computed ther- 
mal pulse. The boxes represent the range of observed 
Se and Kr abundances for the full PN sample, see text 
for details. 



els plotted in Fig. [3] There are still large uncertainties 
in deriving the Kr and Se abundances from PNe spec- 
tra. Hence it is possible that the most extreme PN 
abundances are effected by these uncertainties. 

Finally, one other possibility is that the PN with 
the lowest Se abundances do not reflect internal nucle- 
osynthesis during the progenitor AGB phase, but the 
initial Se abundance. Hence the very low Se/O abun- 
dances are the result of a peculiar galactic chemical 
evolution. If this is the case, then it would be intrigu- 
ing to see if chemical evolution models would be able to 
explain the large spread in Se found in the full sample 
of PNe. 

4 Low-metallicity PN 

There have been a few low-metallicity PN identified 
inclu ding BoBn 1 (|Otsuka et alj 120081 ; IZiilstra et all 
2006). These PN likely originate from low-mass pro- 
genitors with initial masses less than ~ IMq. These 
PN could be used as tracers of the low-metallicity en- 
vironment out of which the progenitor stars formed, 
and they could be used to constrain the evolution and 
nucleosynthesis of low-metallicity AGB stars. 

In Fig. |3 we show the s-process abundance distri- 
bution for two low-metallicity, low-mass AGB mod- 
els. For these models we choose an initial a-enhanced 
abundance pattern (e.g., [O/Fe] = +0.4), which re- 
sults in initial heavy element abundances with [X/O] 
= —0.4 when measured relative to oxygen. Note that 
relative to iron, the initial abundances for these ele- 
ments are scaled solar, i.e., [X/Fe] = 0. In Fig. [4] we 
show the observed Se and Kr abundance distribution 
from the SD08 sample, noting that most of the PN 
used were Galactic PNe with higher metallicities than 
used in the AGB models presented here. From Fig. [3] 



4 



Publications of the Astronomical Society of Australia 



CN - 




30 35 



Proton number, Z 

Figure 4: Same as Fig. [2] except the models shown 
are the L25M , [Fe/H]= -2.3 model (solid line), and 
the 2M , [Fe/H]= -2.3 model (dashed line). An 
a-enhanced initial abundance pattern was used (e.g., 
[O/Fe] = +0.4), and scaled solar for elements heavier 
than iron. 



we see that the relative ratio of Zn/O is sub-solar, with 
[Zn/O] ~ -1 (starting at [Zn/O] = -0.4) from both 
models. If we examine the oxygen abundance, we find 
that a significant amount of this element is dredged 
to the surface, that is, [O/Fe] = 0.87 for the 1.25M 
model, and [0/Fe]= 1.20 for the 2M case. In both 
cases, Zn is actually produced by the s-process but in 
smaller quantities to oxygen, hence the [X/O] ratios 
do not reflect the degree of production. For example, 
the final [Zn/Fe] = 0.45 for both models. Note that in 
both models, the iron abundance remains essentially 
unchanged. 

Likewise for Kr, the abundances measured rela- 
tive to iron are [Kr/Fe] = 1.0 dex and 1.1 dex for 
the 1.25M0 and 2Mq models, respectively. This is 
compared to [Kr/O] ~ (see Fig. [4}. Hence our re- 
sults indicate that for the lowest metallicity PNe, the 
dredge-up of oxygen implies that this element is no 
longer a suitable proxy for the initial metallicity of 
the progen itor star. Argon is somet imes used in place 
of oxygen (jLeisv fc Dennefeldl[2006l . SD08), under the 
assumption th at it remains unchange d by AGB nucle- 
osynthesis. In Ka rakas et al.l (j2009h we showed that 
this is indeed the case. Zinc is also a potentially use- 
ful reference element beca use it does not con dense into 
dust as easily as iron (e.g..lWeltv et alJll999h. and can 
be observed in PNe (jDinerstein fc Geballdl200lh . Zn 
is at the beginning of the s-process chain, and our re- 
sults show that it can be produced in observable quan- 
tities in low-metallicity AGB stars. The best indi- 
cator of the metallicity is iron, but the element iron 
abundance is difficult to accurately determine from 
PN spectra ilPerinotto et al.ll 19991 ; ISterling et aL1l2005l : 
iDekado Inglada et al.ll2009h . 

5 Concluding Remarks 

Abundances derived from PN spectra are an invaluable 
tool to study the late phases of stellar evolution of low 
and intermediate-mass stars, as well as trace the chem- 



ical evolution of galaxies. Elemental abundances for C, 
N, O, S, the noble gases He, Ne, and Ar, as well as el- 
ements produced by neutron capture processes includ- 
ing Ge, Se, Kr, and Ba can be obtained, with varying 
degrees of accuracy. Observations of neutron capture 
elements in PNe can be used to constrain the amount 
of third dredge up mixing follo wing the final therma l 
pulses of the progenitor star. In iKarakas et al.l (|2009f ) 
we attempted to use the Se and Kr abundance informa- 
tion for Type I PNe, along with information on lighter 
element abundances (e.g., C, N, O, Ne etc) to con- 
strain the neutron source operating in intermediate- 
mass AGB stars. We were unable to reach firm con- 
clusions using the Se and Kr abundances. Certainly, it 
seems that the 22 Ne source is required to produce en- 
richments in s-process elements in these more massive 
AGB stars. It is less clear if the 13 C neutron source 
is also required. Increases in the Kr abundances were 
observed in models with a 13 C pocket, but not beyond 
the amounts observed in Type I PN spectra. Obtain- 
ing abundances for more heavy elements beyond the 
first s-process peak may help to distinguish among the 
possibilities. 

Type I PNe may also originate from binary interac- 
tions, and/or by rapid rotation during the progenitor's 
main sequence phase. Certainly obtaining more abun- 
dances for Type I PNe may help to distinguish between 
the single or binary evolutionary channels, and to con- 
strain the initial progenitor masses of these systems. 

We compared results from lower mass AGB models 
to the full sample of PNe from SD08. Our results indi- 
cate that AGB models of solar or near solar metallicity 
produce a reasonable match to Se and Kr enrichments 
seen in the full sample. Our models cannot explain 
the PNe with sub-solar Se abundances, nor the PNe 
with the largest Kr enrichments. These outliers are 
only a small percentage of the full sample (~ 5%), and 
may reflect either uncertainties in the abundance de- 
terminations or the effects of inhomogeneous chemical 
evolution. 

Finally, we present results for low-mass, low metal- 
licity AGB models with metallicities appropriate for 
the Halo ([Fe/H] « -2.3). These AGB models were 
computed with an a-enhanced initial composition, and 
scaled solar for elements heavier than iron. In these 
low-metallicity AGB models, oxygen is dredged-up to 
the surface, which means that measuring the enrich- 
ment of heavy elements relative to this element (e.g., 
[Zn/O]) produces misleading results. For example, 
[Zn/O] » -1 for these models, when [Zn/Fe] « 0.45. 
Hence, it is important to determine the level of O 
dredge-up occurring, or to use another, more suitable 
proxy for the metallicity (e.g., Ar or Fe if possible). 
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